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We identify and study new nonlinear axisymmetric equilibria with incompressible flow of arbi¬ 
trary direction satisfying a generalized Grad Shafranov equation by extending the symmetry analysis 
presented in [G. Gicogna and F. Pegoraro, Phys. Plasmas 22 , 022520 (2015)]. In particular, we con¬ 
struct a typical tokamak D-shaped equilibrium with peaked toroidal current density, monotonically 
varying safety factor and sheared electric field. 

PACS numbers: 


The analysis of symmetry properties of ordinary or 
partial differential equations is a very useful and fruit¬ 
ful tool for studying the general structure and the space 
of solutions and for finding explicit solutions. We con¬ 
sider symmetries described by continuous Lie groups of 
transformations, see for example i-i- In particular 
such symmetry techniques were applied to construct lin¬ 
ear and nonlinear solutions of the Grad-Shafranov equa¬ 
tion i-d. 

This paper is concerned with the study of a class of 
solutions obtained through Lie-group-symmetry analy¬ 
sis of the following generalized Grad-Shafranov equation 
(GGS) @-[ni: 
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Here, the poloidal magnetic flux function '(/)(r, z) labels 
the magnetic surfaces, where (r, </>, z) are cylindrical co¬ 
ordinates with z corresponding to the axis of symmetry; 
Mp{ip) is the Mach function of the poloidal fluid velocity 
with respect to the poloidal Alfven velocity; X{'tp) re¬ 
lates to the toroidal magnetic field, = I/R, through 
I = X/{1 — M^); $(■!/') is the electrostatic potential; 
for vanishing flow the surface function Ps{ip) coincides 
with the pressure; B is the magnetic field modulus which 
can be expressed in terms of surface functions and R', 
A* = • iy/R^)] and the prime denotes derivatives 

with respect to ijj. Because of incompressibility the den¬ 
sity p{ip) is also a surface quantity and the Bernoulli 
equation for the pressure decouples from (0: 


P = P,{^)-p 
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where v is the velocity modulus. The quantities Mp{ip), 
X{ip), PsW, pitp) and $(^) are free functions. Deriva¬ 
tion of m and dSDis provided in i-im. 


Eq. ([T]) can be simplified by the transformation 
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under which © becomes 
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Note that no quadratic term as |Vitp appears anymore 
in (III). It is emphasized that once a solution of (j4]) is 
obtained, the equilibrium can be completely constructed 
with calculations in the u-space by employing ([3]), and 
the inverse transformation 

pU 

^{u)= [l-Ml{f)]-^^^df (5) 
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For example, one has for the electric field 


E = -V$ = - 


dip 


du 


Xu 


Contrary to what is stated in [T^ the (explicit) inversion 
Ip = if {u) is not needed provided that and the other 
free surface quantities are assigned as function of u. For 
parallel flows ($' = 0), Eq. (jH) reduces in form to the 
usual GS equation. 

The symmetry properties of (HJ, expressed in terms of 
alternative surface quantities (Eqs. (13)-(16) of l^) in 
connection with the variational derivation of [l4|, were 
studied in [T^ (see also @ and [il). Henceforth, we 
follow closely the analysis developed in [I^ and extend 
the results therein. 

Choosing the free function terms in (j4]) as 
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where oq, 02 , 04 are free parameters, 0 assumes the form 


(7) 
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This equation admits the Lie point scaling symmetry 
^ d d 1 d 

and the ‘exceptional’ symmetry 
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If y = — crr^ is introduced as a weaker type ol 

conditional symmetry (see [^, [l^ l then we can write 
the GGS equation in terms of the T—invariant variable 
s = ar^ + (cr 7 ^ 0,1). Then the above mentioned 
symmetries map solutions into solutions of the form 


y\ ■= 


r 

SA 


u{r, z) = ag, 02 , ^ 4 ) 

sa := crr^ + [z + A(r^ + z^)]^ 
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For cr = — 1, (TO recovers the class of solutions ob¬ 
tained in [T^ without the restriction of constant den¬ 
sities adopted therein; solutions of the form (ITU)) hold foi 
arbitrary Mach functions, Mp(it) and densities p{u). 

We will construct and study one tokamak pertinent 
solution of this symmetry-generated class. Substituting 
into © we obtain 


[4cr(cr - l)y^-I-4(1 - 2cr)y^-I--|- 
-|-[6cr(cr - 1 ) 2 /^ -I- (4 - 6cr)y](j) + 
1 3cr(l-cr) ao a^y , a^y"^ 
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In order that the coefficient of (f) in vanishes for 
a particular y = yo, we choose yg = lj{cr — 1). Then, to 
integrate m we employ the initial conditions ipiuo) = 1 
and 


Tyo*/* (yo) = 


, 1 -I- 3ct 


-I- Oq + Cl2yo + Cl4yo] (12) 


the latter one stemming from (HU for y = yg. We solve 
numerically (ED by choosing cr = 1/4, A = a/3/4, ag = 
1/4, a 2 = —0.5, 04 = 0.1 and integrating in the interval 
yo ^ y ^ yoj with yg = 1/cr. Using cubic fitting we found 
that 


^(j;) = -1.341?/3 - 4.343y2 - 4.999y - 1.106 (13) 

Thus, we have obtained the equilibrium configuration 
with a crescent-shaped cross-section shown in Fig. 1. 
By exploiting the invariance of GGS under constant dis¬ 
placements z —>■ z -I- d, the configuration has been prop¬ 
erly displaced along the z-axis so that the magnetic axis 
lies on the plane z = 0. The outermost magnetic sur¬ 
face corresponding to the flux value of Uout = —125 Wb 


touches the z-axis at a couple of corners. At the magnetic 
axis located at rg = 1.19, Zg = 0.05 we have ug = —135 
Wb. A more pertinent tokamak equilibrium can be con¬ 
structed by choosing a fixed boundary to coincide with 
an interior D-shaped magnetic surface. Such a boundary 
corresponding to the value of Ub = —132.543 Wb is indi¬ 
cated by the black-colored curve in Fig. 1 and separately 
in Fig. 2. 



(r-1.2) 


FIG. 1: u-contours on the poloidal plane for an equilibrium 
of the class of solutions |IQ|. The black-colored D-shaped 
curve corresponds to the boundary of a tokamak pertinent 
equilibrium shown separately in Fig. 2. 

To completely construct the equilibrium we choose the 
Mach-function and density peaked on the magnetic axis 
and vanishing on the boundary as M‘^{u) = Mq (1 — 
M^/uj), p = po(l — l^b) respectively, where 
and po = 4 X 10“^ Kg/m^. The pressure of the D-shaped 
equilibrium of Fig. 2, peaked on axis and vanishes on 
the boundary, is shown in Fig. 3. The safety factor has a 
typical tokamak variation monotonically increasing from 
the magnetic axis to the boundary, as can be seen in Fig. 
4. The toroidal current density is also peaked on axis and 
vanishing on the boundary (Fig. 5) unlike a hollow 
profile of an equilibrium with parallel flow correspond- 
ing to the same ansatz with 04 = 0 constructed in 
pU| . This result indicates that the electric field associ¬ 
ated with the non-parallel component of the flow, shown 
in Fig. 6 , may play a role in determining the equilibrium 
characteristics. 

In summary, by employing Lie and weak conditional 
symmetries of a GGS equation for plasmas with incom¬ 
pressible flow of arbitrary direction we have extended a 
class of nonlinear axisymmetric solutions obtained in [iU| . 
In particular, we constructed a D-shaped configuration 
with sheared electric field and typical tokamak equilib¬ 
rium characteristics, i.e. pressure and toroidal current 
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FIG. 2: M-contours on the poloidal plane for the tokamak 
D-shaped equilibrium. 



FIG. 3: The pressure profile on the plane 2 = 0 for the D- 
shaped equilibrium of Fig. 2. It vanishes on the boundary 
of the D-shaped equilibrium at Ub = —132.543 Wb. Here 
Pmax = 0.316 Atm. 


density peaked on the magnetic axis and safety factor 
monotonically increasing from the axis to the boundary. 
It would be interesting to extend the search for other 
relevant equilibria by potentially identifying other weak 
conditional symmetries following the procedure of @ and 
[T^ . Work along these lines is in progress. 



FIG. 4: The safety factor for the D-shaped equilibrium of Fig. 

2 . 



FIG. 5: The current density for the D-shaped equilibrium of 
Fig. 2. It vanishes on the boundary of the D-shaped equilib¬ 
rium. Here Jmax — l-TMAjm^. 
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FIG. 6: The electric field for the D-shaped equilibrium of Fig. 
2. It vanishes on axis. Here Emax = 1 Kv/m. 
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